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Anderson localization is a paradigmatic coherence effect in disordered systems, often analyzed 
in the absence of dissipation. Here we consider the case of coherent dissipation, occuring for open 
system with coupling to a common decay channel. This dissipation induces cooperative Dicke super- 
and subradiance and an effective long range coupling, expected to destroy Anderson localization. We 
are thus in presence of two competing effects, i.e localization driven by disorder and derealization 
driven by dissipative opening. Here we show that in an open 3D Anderson model, subradiance 
enables the system to preserve signatures of a localization transition. We demonstrate the existence 
of a subradiant localized regime, emerging from the interplay of opening and disorder, in which 
subradiant states are hybrid and localized, while superradiant states are extended. We also provide 
analytical predictions for this regime, confirmed by numerical simulations. 

PACS numbers: 72.15.Rn, 03.65.Aa, 03.65.Yz, 05.60.Gg 



I 

in 



c3 



1 

C 

o 
o 



in 
in 

o 



X 

S-H 



One of the most important effect induced by quan- 
tum coherence is Anderson localization [JJ. This effect 
is relevant for many areas of Physics, quantum comput- 
ing [21 [3], cold atoms [4], light harvesting systems [5], 
mesoscopic electronic devices [BJ. Initial search of ex- 
perimental evidence of Anderson localization for non in- 
teracting waves in 3 dimensions has been limited by the 
presence of absorption [7j. Recent progress has how- 
ever been obtained in a large variety of fields of re- 
search, including acustic waves [5], light waves [5], matter 
waves [TU] . 

Whereas absorption and dissipation is often considered 
to be limiting phenomena for observing Anderson local- 
ization, the situation of opening with corresponding co- 
herent dissipation has not been systematically addressed. 
This situation has been considered in a different research 
community [SJ [IT], investing cooperative effect induced 
by coupling to a commun decay channel, following the 
pioneering work by Dicke in 1954 [12]. The question 
whether localization can survive coherent dissipation in 
open 3-D system has been discussed in particular in the 
case of resonant light scattering [U [T3] . 

In this letter we focus on the case of coherent dissi- 
pation due to the coupling to a common decay channel, 
relevant for many realistic situations [jTHHUJ] ■ If this cou- 
pling is strong enough, this kind of dissipation induces a 
well known important coherence effect: Dicke super- and 
subradiance [12] . which induces the so-called Superradi- 
ance Transition (ST), driven by the opening, i.e. by the 
coupling to an external environment characterized by a 
continuum of states. This transition is to be compared 
on the other hand to the Anderson localization, driven 
by intrinsic local disorder which consists in the suppres- 
sion of diffusion due to exponential localization of all the 
eigenfunctions of the system. 

Both Dicke superradiance and Anderson localization 



have been widely studied in literature in a separate way, 
and their interplay has been poorly analyzed. In order 
to highlight the fundamental aspects of this interplay, 
we have choosen to analyze a simple but general model: 
the 31? Anderson model [JJ. A related ID model, with- 
out a metal-insulator transition, has been already con- 
sidered in Rcf. 15j. A common feature of disordered sys- 
tem with coherent dissipation is the competition between 
long range hopping induced by opening and localization 
induced by disorder. It is commonly accepted that any 
kind of long range hopping, decaying slower than l/r d , 
where d is the system dimension, destroys localization [lj. 
Here we show on contrary that the very correlated na- 
ture of long range hopping due to the opening induces a 
Superradiant Transition which allows to preserve local- 
ization. 

Dicke superradiance occurs in the large opening regime 
and it cannot be treated by perturbation theory. In such 
a case the effective non-Hermitian Hamiltonian approach 
to open quantum systems has been shown to be very 
effective [TB]. Non-Hermitian Hamiltonians have been 
already employed in random matrix theory |17H20) . in 
paradigmatic models of coherent quantum transport |21[ 
122"] and in realistic open quantum systems [U UU r2"5H2"5] . 

In the closed 3D Anderson model, a particle hops be- 
tween neighbors sites of a 3D cubic lattice with N sites, 
in the presence of on-site disorder. The Hamiltonian of 
the closed 3D Anderson model can be written as: 



N 



ff = ^£; j |i)0'| + n^(|i)(i| + |i)oi) 



(i) 



J'=l (i,j) 



where the summation (i,j) 



runs over the nearest- 



neighbor sites, Ej are random variables uniformly dis- 
tributed in [— W/2,+W/2], W is a disorder parameter, 
and 51 is the tunneling transition amplitude. The nature 
of the eigenstates of the 3D Anderson model depends on 



the degree of disorder: for small disorder the states in 
the middle of the energy band are extended, while close 
to the band edges, for energies below the mobility edges, 
they are localized [25] • On increasing W/Cl, the mobility 
edges approach one to each other and above the critical 
value W/Q, ~ 16.5 26 , all states became localized and 
a global AT occurs. In the localized regime the shape 
of eigenfunctions behaves as \ip(j)\ ~ exp(— \xj — £o|/£), 
where Xj are position vectors and £ is the localization 
length. In this work we will disregard the effect of the 
mobility edges, focusing on the global AT: 



AT at W/Q = 16.5 



(2) 



Our main question is about the effect of opening on 
the localization transition which occurs at AT. We open 
the 31? Anderson model by allowing the particle to es- 
cape the system from any site into the same continuum 
channel. This situation of "coherent dissipation" , can be 
met in many realistic systems [3J H] , when the wavelength 
of the particle in the continuum channel is comparable 
with the sample size. The case of only one channel in the 
continuum is somehow extreme, but we expect that our 
findings have general validity whenever more states com- 
pete to decay in the same channel. This kind of opening 
is different from a standard coupling to a thermal bath 
(for instance described by time dependent diagonal terms 
in the Hamiltonian [57]), and takes into account only 
the particle escape from the system (dissipation). The 
open system is described by the effective non-Hermitian 
Hamiltonian [T3] : 



(H t 



eSjij 



(Ho); 



■1 , 



i-ji;w=(%-^. 



(3) 



where Hq is the Anderson Hamiltonian, Eq. (fTl), A\ are 
the transition amplitudes from the discrete states i to 
the continuum channels c. In our case we have one decay 
channel, c — 1, equal couplings ( A\ = ^/j) so that 
Qij = 1 Vi,j is a full matrix. The complex eigenvalues of 
-ffeff can be written as, E r — 
widths of the states. Since the average width is given 
by 7 |16j , the degree of resonance overlap is determined 
by k = j/D, where D is the mean level spacing of the 
energy levels of Hq. We can regard k as an effective 
degree of opening and at k ~ 1 [El EH] a segregation 
occurs, i.e. almost the entire decay width is allocated to 
just one short-lived "superradiant" state, while all other 
N — 1 states become "subradiant" with a small decay 
width. We refer to this segregation as the ST. We note 
that our k parameter bears some resemblance to a often 
used Thouless parameter g, which in the case of opening 
at the edges of the system has been shown to provide 
relevant information on a metal- insulator transition |28j . 
Here, the opening is obtained by coupling all sites to a 
commun decay channel, and our n thus does not reveal 
the sensitivity to the boundaries of the system. 

The presence of disorder affects ST since D depends 
on disorder; indeed for small disorder, the width of the 
energy band is given by $7 and D ~ 1251 /N [JJ, while for 
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FIG. 1: Log(PR — 1), where PR is the participation ratio, 
is plotted in the 7/fi and W/Q plane. As a dashed line AT 
is shown, Eq.|2|, while the full line refer to the ST at k — 1, 
Eq.pl), so that the region above it re > 1 is characterized 
by superradiant and subradiant states, differently from the 
region below it (re < 1). For an explanation of the different 
regimes see text. Panel (a) refers to the state with the largest 
width, which become superradiant above the ST. Panel (b) 
refers to the other N—l states which become subradiant above 
the ST. In both panels the system is a cubic lattice with TV = 
10 x 10 x 10 sites and £2 = 1. Each data is obtained averaging 
over 5 realizations of disorder, for panel (b) an additional 
average over all states has been done. 



where T r are the decay large disorder, D ~ W/N, so that, 



ST at K : 



2L 
D 



1 with 



-7JV/12Q if W/n < 12 
-77V/W ifW/n»12 

(4) 
So, above ST (k > 1), we have a superradiant regime, 
with superradiant and subradiant states, while below 
ST(k < 1), all states are affected by the opening in a 
similar way. From Eq.Q it is clear that we can control 
the effective degree of opening k by varying the strength 
of the disorder. 

Since Q is a full matrix the opening induces an effec- 
tive long range hopping among the sites, which is gen- 
erally expected to destroy localization [35]. The long 
range hopping can be explained as a second order effect 
of coupling via a commun vacuum mode which results in 
the fact that the particle escaping from one site can be 
reabsorbed in far sites before leaving the system. This 
interaction mediated by the continuum, is at the origin 



of the ST. Thus disorder and opening have opposing ef- 
fects: while disorder tends to localize the eigenfunctions, 
the opening tends to delocalize them due to the induced 
long range hopping. 

In order to analyze the interplay of disorder and 
opening we first study the participation ratio, PR = 
(Et K^)! 4 ]" 1 ) of the eigenstates \ip) of H cS in Eq. J3J), 
where (. . .) stands for the average over disorder. The 
PR is widely used to characterize localization proper- 
ties [30] : for extended states it increases proportionally 
to the system size, while it is independent for localized 
states. Note that since PR > 1 in the following we will 
analyze PR — 1. 

This interplay between disorder and opening generates 
a rich phase diagram, as shown in Fig. [T] where we plot 
Log (PR — 1) in the 'y/Q-W/fl plane for a system with 
N = 10 3 sites. In the same figure we also show AT 
(dashed curve) and ST (k — 1 full curve). In the upper 
panel of Fig. [T] we analyze the state with the largest 
width which become superradiant above ST, while in the 
lower panel we considered the other N — 1 states, which 
become subradiant above ST. As one can see the response 
to disorder in the superradiant regime is very different for 
superradiant and subradiant states: while the former do 
not feel AT and remains delocalized up to ST, the latter 
are sensitive to AT. The reason of this different behavior 
can be explained from perturbation theory and it is due 
to the fact that the superradiant state is very far from 
the complex energies of the other states and thus it is less 
affected by the disorder [T5] . Below the ST, all states feel 
the disorder and the opening in a similar way. From Fig. 
[I] it is also clear that the parameters which determine 
the nature of the eigenstates are W/VL and the effective 
coupling strength to the continuum k, Eq. Q . From Fig. 
[llwe can distinguish three main regimes: 
Region I (at the left of the AT line in Fig. [I]) : most of 
the states of the closed system are delocalized, and the 
opening does not change their extended nature. For this 
reason we do not distinguish between k < 1 and k > 1 . 
We neglect the discussion about the effect of the opening 
on the mobility edges of the closed system, since we are 
interested in the global transition to localization at AT. 
Note that for large opening and small disorder we recover 
the usual Dicke regime. 

Region II (right of the AT line and below the ST line in 
Fig. II|): the states of the closed system are localized and 
even if the effective opening is small (we do not have su- 
perradiance) , it induces hybrid states (see discussion be- 
low). It is important to note that this region disappears 
for large N values, since, from Eq. Q the critical ST line 
is, 7/O = 12/JV for W/n < 12 and 7/Q = W/(Q,N) for 
W/fl ^12. Note that in the limit of large disorder and 
small opening we recover the Anderson localized regime. 
Region III (right of the AT line and above the ST line 
in Fig. [I]): the states of the closed system are localized, 
but the effective opening is large and we have superradi- 
ance. This regime is very interesting since the superradi- 
ant state is fully delocalized (Fig. [II upper panel), while 
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FIG. 2: 3D Anderson model. Average participation ratio 
(PR - 1} vs N (up to Nmax = 27 3 ). Here 7/fi = 1 and 
W/Q. = 80, so that we are in the Subradiant localized regime. 
Full circles stand for the superradiant state, and open red cir- 
cles stand for subradiant states. Each point has been obtained 
by averaging both over the states and disorder. Dashed lines 
are plotted to guide the eye. 



the subradiant states are hybrid and localized (see dis- 
cussion below and Fig. [3]): The existence of this regime, 
which we call the subradiant localized regime, is the main 
result of this work. 

A full analysis of the three regions will be presented in 
a future work. Here we focus on the subradiant localized 
regime, which is the regime with very large opening, in 
which one might expect no signature of localization. 

In order to establish the localized or extended nature 
of the states in the different regimes, it is necessary to 
analyze the scaling of the PR with the system size. In 
Fig. [2] PR — 1 is shown vs N for the subradiant and the 
superradiant states in the subradiant localized regime, at 
fixed strength of disorder. As one can see the superra- 
diant states (full circles) are extended, PR oc N, while, 
surprisingly, we find that the subradiant states behave 
like localized states with a PR independent of N (open 
red circles). We have also checked the behavior of the 
PR in other regions (not shown): Below AT (Region I) 
we find that the states are extended (PR oc N) , as in the 
closed model. In (Region II), we have found PR oc 1 + N 
which would indicated extended states. In reality the 
limit N — » 00 can not be performed since at N ~ W/j 
we cross the ST, entering the subradiant localized regime 
and the Region //disappears. We can conclude that this 
is a finite-sized effect. 

The behavior of the PR is not enough to character- 
ize the nature of the eigenstates tp, it is also important 
to analyze their structure, by means of the average den- 
sity profile, and defined as (|?A| 2 ). In the subradiant lo- 
calized regime the superradiant states are fully extended 
and are well approximated by the following expression: 
I Superradiant) fa -j^Eil*)) wnere I*) are the site states 
of the Anderson model. The analysis of the subradi- 
ant states in the subradiant localized regime shows that 
they are hybrid states [13] with an exponentially local- 
ized peak plus a constant plateau. The exponentially 
localized peak behaves like a localized state of the closed 
Anderson model, compare respectively solid and dashed 
curves in Fig. [3] On the other hand, the height of the 




10 

m J 

io 2 
£io° 
v io- 2 

io" 4 

10" 



ijOoooooooooooooooooooodoooooooooooooaoo 



Region I 



Subradiant Localized 



AT 



ST 



Region II 



10 

w/n 



FIG. 4: Average width vs W/Q for N = IO 3 , and 7 = f2 = 
1. (Black) circles stand for the width of the superradiant 
state, while (red) crosses stand for the average widths of the 
TV — 1 subradiant states. The vertical dashed line represents 
the AT, while the full vertical line indicates the superradiant 
transition (ST). 



FIG. 3: All data refer to the subradiant localized regime. 
Full curves are for the open system, while the dashed curves 
are for the closed one. In the two upper panels, the average 
density profile of subradiant states as a function of the dis- 
tance r from the heighest Anderson peak is shown. Left upper 
panel: N = IO 3 , 7/fi = 100 and different value of disorder 
as indicated in the legend. Right upper panel: W/Q = 80, 
7/fi = 1 and different system sizes as indicated in the legend. 
Lower panel: A single hybrid subradiant state is compared 
with perturbation theory, Eq. S. 



constant plateau is independent of both W and 7, see 
Fig. [3] (left panel) , and it decreases with the system size 
as 1/JV, Fig. [3] (right panel). 

For large disorder, W 3> f2 and above ST (subradiant 
localized regime) , following Ref. |15L 120) it is possible to 
obtain an analytical expression for the localized subradi- 
ant states: 



Im> = 




with 



JV 



1 



E„ 



= 0. 



(5) 
where C M is a normalization factor, |j°), Ejo are 
the eigenstates and eigenvalues of the closed Anderson 
model, Eq.pT), and e M are defined by the constraint 
given in the second equation of Eq.Q. Note that /i = 
1, .. , JV — 1 spans the subradiant subspace and each e M 
lies between two neighbor levels Ejo . Even if Eq. (l5|) de- 
scribes very well the localized subradiant states, see Fig. 
[3] lower panel, it is not trivial to derive a close expression 
for the density profile. 

We have also computed the density profile of the states 
for small opening and large disorder Region II. When the 
effective opening is small (below ST) , the problem can be 
treated by standard perturbation theory and similarly to 
what was found in |15j . we find that also in this case we 
have hybrid states, with a plateau height proportional 
to (j/W) 2 but independent of JV. The different scaling 
of the plateau height with the system size, is at the ori- 
gin of the different scaling of the PR of the open system 
between Regions II and the subradiant localized regime. 
Indeed for very large disorder, the hybrid states in both 
regions, are mainly localized on only one site with prob- 



ability a and on all the other JV — 1 sites with proba- 
bility 6, and we have £\ |-0i| 4 = a 2 + b 2 (N - 1), with 
a + b(N — 1) = 1. So that when b ex 1/JV (subradiant lo- 
calized regime) , the PR is independent of JV in the large 
JV limit, while, when b is independent of JV (Region II) 
one has PR ~ 1 + 26JV in the limit bN <C 1, in agreement 
with the discussion given above. 

A signature of the AT can also be found in the behavior 
of the decay widths of the subradiant states. In Fig. [4] 
we show the decay width of both the superradiant and 
subradiant states as a function of the disorder. While the 
former is not affected by disorder up to the ST, the latter 
feels the disorder above the AT (note that, increasing 
the disorder at fixed 7, the effective degree of opening is 
decreased and the three regions are crossed). For very 
large disorder (Region II), all widths approach the same 
value To — 7, which corresponds to the decay width of 
an isolated site. 

Conclusions We analyzed coherent effects induced by 
the interplay of coherent dissipation and disorder. For 
this purpose we considered a 3D open Anderson model 
in which a particle can escape from any site to a com- 
mon channel in the continuum. This kind of opening 
induces a strong long range hopping (all-to-all) between 
the sites of the Anderson model. Contrary to expec- 
tations, we show that the opening does not destroy all 
features of Anderson localization. Remarkably, for large 
effective opening, where we have superradiance, we estab- 
lished the existence of a subradiant localized regime with 
extended superradiant states, and localized subradiant 
states (with a size independent participation ratio). We 
determined, both numerically and analytically, that the 
subradiant localized states have an hybrid nature, with 
an exponentially localized peak and a constant plateau. 

This surprising result can be only explained by the 
highly correlated nature of the long range hopping 
present in this kind of systems. These correlations induce 
a Superradiant transition where the coupling with the ex- 
ternal world is taken by the superradiant state, leaving 
the subradiant states effectively decoupled, and thus able 
to preserve the localized nature of the closed system. Let 
us observe that recently in Ref. [3T] , localized states have 
been found in a related model with hermitian long range 



hopping, thus supporting the generality of our results. 

Our analysis is relevant both from an applicative point 
of view, in the search of Anderson localization in 3-D dis- 
sipative systems (i.e. cold atoms), and from a theoretical 
point of view (role of long range hopping in the metal- 



insulator transition). 
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